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Abstract We present results of models of the physical 
space and parameters of the accretion disk of Sagit- 
tarius A*, as well as simulations of its emergent spec- 
trum. This begins with HARM, a 2D general relativis- 
tic magneto-hydrodynamic (GRMHD) model, specifi- 
cally set up to evolve the space around a black hole. 
Data from HARM are then fed into a 2D Monte-Carlo 
(MC) code which generates and tracks emitted photons, 
allowing for absorption and scattering before they es- 
cape the volume. 
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1 Introduction 

Sagittarius A* is the supermassive black hole (about 
3.6 million Solar masses) at the center of our galaxy, 
which is known as a strong radio source with specific 
x-ray signatures which have come to be associated with 
flaring or quiescent states. The short time scales of 
these flares, as well as the lack of significant time lags 
between NIR and x-ray components, suggests the pos- 
sibility of these arising in a small volume close to the 
event horizon. 

Recent attempts to model the spectrum of Sagittar- 
ius A* have proven promising, but have relied upon 



simplifications or assumptions which have somewhat di- 
minished the impact of these results. Our intent is to 
deliver the most consistent modeling possible with cur- 
rent technology, and avoid some of the simplifications 
that previous simulations have required. 

This begins with using a GRMHD code specifically 
designed to evolve an accretion disk around a black 
hole. This is a considerable advantage over models 
which specify analytic formulae to describe density, 
temperature, and magnetic field profiles. Even simu- 
lations done with MHD codes suffer some of the same 
downfalls as simpler set-ups — specifically that effects 
of general relativity near a supermassive black hole be- 
come quite important. 

Data output by the GRMHD code are used as input 
for our MC code. Unlike other models which rely only 
on emission and scattering calculated analytically, this 
simulation creates a realistic situation through use of 
the Monte-Carlo random generation method for pho- 
tons, and allows for scattering and absorption before 
photons escape the volume. 

It should be noted that both of the codes used here 
are 2D and axially symmetric. The third dimension 
would be useful for full consistency, but the authors feel 
that its lack of inclusion does not preclude the results 
from being significant. 
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2 HARM GRMHD Code 

Our GRMHD modeling is done using the HARM (High 
Accuracy Relativistic Magnetohydrodynamics) code 
developed by Gammie et al. ( 2003 1 . The code simulates 
GRMHD evolution through time by the use of particle 
number conservation, the four energy-momentum equa- 
tions, the MHD stress-energy tensor, and the induction 
equation. These equations take conserved hyberbolic 
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forms, for easy integration: 
dF{U) 



dU 
~dt 



^S{U) 



(1) 



where U, F, and S represent the set of conserved vari- 
ables, fluxes, and sources, respectively. 

HARM is set to first seed a torus around a black hole 
with an initial density and poloidal magnetic field, then 
evolve the space through time, tracking density, inter- 
nal energy, 4-velocity, 4-magnetic field. This is allowed 
to continue until the space reaches a "steady state" - of 
special note in the code output is the highly turbulent 
magnetic field, an artifact of the magneto-rotational in- 
stability, which is an important mechanism governing 
the transport of angular momentum in the disk. This 
instability causes the turbulence which is likely to be at 
least partly responsible for the energizing of electrons 
in the disk, though in our MC trials, we only consider 
fully thermal distributions. 



3 MC Emission and Scattering Code 

The MC code used for photon emission and scattering 
has a long history and is discussed in a number of re- 



sources, such as 


Canfield et al. 


(1987 


); Liang & Dermer 


(1988 


1; 


Bottcher et al. 


(1998) 


For a complete treat- 



ment, readers should see these papers. 

In general, this code is a coupled MC/FP (Fokker- 
Planck) code. For our intents, the FP evolution of the 
electron distribution was unnecessary at this stage, so 
it was turned off to allow a fixed temperature given by 
the HARM output. The code is set up on a 2D axially- 
symmetric cylindrical grid, creating a (hollow or solid, 
depending on whether the inner radius is set to zero) 
cylindrical shape. Each zone is assigned a density, ion 
and electron temperatures, magnetic field amplitude, 
and thermal and nonthermal distribution components. 
For our purposes, this is simply set to be a Maxwellian, 
but the code allows power law nonthermal distributions 
as well. The code allows emission from the volume and 
boundaries, and emitted photons are tracked and al- 
lowed to scatter or absorb. 

Output from the HARM code is used to assign the 
zone quantities. That is, we specify a maximum den- 
sity and the rest of the output scales to give us the ion 
temperature and saturated MRI magnetic field. As- 
suming that inclusion of an initial toroidal magnetic 
fields would add to the final field, as this dimension 
is mostly unaffected by evolution, allows us to set the 
magnetic field, as long as it doesn't drop below the 
saturation value. Similarly, as we have a specified ion 
temperature but no way to directly evaluate the elec- 
trons' acceleration, we set a global ratio between the 



two values. That is, in one trial the electron tempera- 
ture in a zone could be set to always be twice the value 
of the ion temperature in the same zone. 

3.1 Flaring Results 

For all of the fits presented in figures below, the open 
circles are data points, triangles are upper limits, and 
the bowties denote the Chandra-obtained fiaring and 
quiescent x-ray data points and slopes. 

The first three figures show fits to the flaring x-ray 
point. As shown, the flaring bowtie has almost zero 
slope associated with it. In order to match this, there 
were several options. Obviously, the bremsstrahlung 
component is very flat out to the high energy cut-off 
point. It is also possible for Compton-scattered com- 
ponents to reach these energies and have flat spectra 
- trials are shown for flts using first the first-scattered 
bump, and secondly, the bump arising from photons 
that scatter twice off of hot electrons. 

These three spectral components provide nearly 
equally-adequate fits, if only considering the flaring be- 
havior, though the second Compton trial in Figure 3 
seems to best flt the lower energy components. 

It should be noted that effort was not focused on 
fitting the very low energy radio points. It is expected 
that these arise primarily from a very large volume of 
low density, temperature, and magnetic field around 
the accretion disc, which is beyond the scale of this 
experiment, as our volume cuts off at r = 40M. 

3.2 Quiescent Results 

Once the flaring trials were completed, consideration 
had to be given to fitting the quiescent points with a 
consistent method, given the flaring fits. It was found 
that very good fits to the quiescent point were possible 
by simply dropping the density or temperature inde- 
pendently, for the two Compton bump trials. Unfortu- 
nately, the bremsstrahlung trial did not flt the quiescent 
point regardless of how its parameters were modified, 
as it will always have a very flat spectrum at the x-ray 
point. 

The second Compton bump quiescent flt was achieved! 

by changing only the density in the trial. This is repre- 
sentative of a global mass accretion rate change, which 
could certainly be possible for observed variability from 
Sgr A*. 

Conversely, it was found that for the flrst Comp- 
ton bump trial, dropping the temperature allowed for 
a good fit. This could indicate that the flaring state 
arises due to some large-scale magnetic energy dissipa- 
tion event, such as rapid reconnection over a large area. 
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Fig. 4 Fit to quiescent data with second Compton bump 
component. Maximum values of the 95x95 ceU grid are 
scaled to n = 6.9x10'^ particles/cm^*, T = 8.2x10^ keV, B = 
1.6x10^ G. Total luminosity is 8.2x10^^ erg/s. 
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Fig. 2 Fit to flaring data with second Compton bump com- 
ponent. Maximum values of the 95x95 cell grid are scaled to 
n = 3.4x10* particles/cm^, T = 8.2x10^ keV, B = 1.6x10^ 
G. Total luminosity is 7.5x10^" erg/s. 
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Fig. 5 Fit to quiescent data with first Compton bump 
component. Maximum values of the 95x95 cell grid are 
scaled to n = 5.2x10*' particles/cm^, T = 4.7x10** keV, B = 
13 G. Total luminosity is 3.4x10^'' erg/s. 



Fig. 3 Fit to flaring data with first Compton bump com- 
ponent. Maximum values of the 95x95 cell grid are scaled 
to n = 5.2x10*^ particles/cm^ T = l.SxlO'^ keV, B = 13 G. 
Total luminosity is 2.6x10^'^ erg/s. 
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3.3 MC Code Modifications 

Since the initial project on fitting data from Sagittar- 
ius A* was completed, work has begun on modifying 
the Monte Carlo code to allow it to better undertake a 
number of problems. 

First in this endeavor was rewriting the synchrotron 
emission routines of the code to support the inclusion 
of anisotropic magnetic fields. As one may expect, in 
supermassive black hole accretion discs, the relative 
strengths of the components of the field differ greatly 
- being heavily biased toward the toroidal component. 
Therefore, it is inconsistent to consider the emission 
from synchrotron radiation to be isotropic. It should, 
instead, be heavily angle dependent. Tests with this 
addition showed it to work as intended, and future 
projects will all include this feature. 

As these simulations deal with situations in which 
individual regions of the volume may be moving with 
significantly different bulk velocities, work has begun on 
the inclusion into the MC code of diffcrcmtially moving 
volume elements. This will allow for another degree 
of consistency in results, by appropriately considering 
the boosting and Doppler shifting of photons emitted 
in each region. 

4 Conclusions 

Fits to the broadband spectra of Sgr A* have been 
shown using a combination of GRMHD accretion disc 
evolution and MC radiation transport codes. Two sets 
of these results - those for the first and second Comp- 
ton bump fits - show good fits to the quiescent x-ray 
point, and both sets, with reasonable changes in the 
physical parameters, fit the flaring x-ray point. These 
changes arc indicative of the near-global events, such 
as mass accretion increases, or the rapid dissipation of 
magnetic energy, which would be necessary to explain 
the observed variability. 

The second Compton bump trials are presented as 
the most promising and complete detailed here. The 
fact that these trials can closely match the IR and op- 
tical data, while fitting the x-ray points, makes these 
the most consistent with observations. 
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